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PREFACE 


This final report under NASA NSG 3109 ME grant is composed of two parts. 

Part I discusses theory relative to measuring inplane displacements 
using speckle interferometry. From the inplane displacement field, 
strains and stresses can be determined. This approach is potentially 
preferred over holographic techniques for evaluating stresses in 
rotating turbine or compressor blades because of the relative ease in 
collecting data, because the rigid body motion requirements are relaxed 
and because data interpretation is simpler. 

Part II discusses theory relative to measuring derivatives of dis- 
placement fields by shearing speckle interferometry. The derivatives 
of displacements are closer to and in some cases are the strain fields. 
Therefore, stresses can be calculated directly. Again this approach *s 
potentially preferred over holographic techniques as noted above. 
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I. INTRODUCTION 

Conventional methods of measuring surface strains and displace- 
ments utilize strain gages, dial gages and various other mechanical 
and electrical sensing devices. The major limitation to these systems 
is the information obtained is for a limited region only. In order to 
determine the displacement or strain field for the entire object, a 
large number of separate measurements must be obtained which is gener- 
ally time consuming, costly and, in some instances, critical areas may 
be missed. 

In recent years, optical techniques employing a coherent light 
source (laser) for detecting and measuring the components of surface 
displacement and strain have come into use. The primary advantage of 
these techniques is that the entire displacement or strain field can be 
determined completely. The majority of these techniques fall into 
three categories: Moire' gauging, holography and laser speckle inter- 

ferometry. The method employed in this report is laser speckle inter- 
ferometry. 

Moire' technique consists of attaching or projecting a grid 
pattern on an object, then photographing the surface of the object 
before and after deformation, and observing the interference fringes 
produced by the overlapping grid pattern [1, 2]*. This particular 
method has the disadvantage that, somehow, a fine grid pattern must 

* Numbers in the brackets refer to the references listed at the end. 
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be generated and superimposed on the test surface and the photographic 
system used to record the image must be of a high sensitivity to be 
capable of resolving the grid system. 

Ennos [3] measured in-plane surface strain using holographic inter- 
ferometry. The principal advantage of this method is the fact that it 
is a very sensitive method of measurement, and the whole surface of the 
structure can be investigated at the same time, rather than a podnt -by- 
point. Although well suited for measuring normal movements, there is 
no general way of eliminating the effect of normal movements from in- 
plane movements. Also, holographic techniques have the disadvantage 
that they require several separate views of the holographic fringe 
patterns of the surface, which requires a considerable amount of data 
reduction to separate out the in -plane displacement field. 

Laser speckle interferometry is a relatively new experimental 
technique which shows promise of alleviating many difficult problems 
in experimental mechanics. The method utilizes simple high-resolution 
photographs of the surface which is illuminated by coherent light. 

The result is a real-time or permanently stored whole-field record 
of interference fringes which yields a map of displacements in the 
object. Suggestions for the direct use of coherent light, in displace- 
ment metrology and contour mapping first appeared in 1968 [4]. Burch 
and Tokarski [5] showed that if two identical speckle patterns are 
superimposed on a photographic plate translated laterally by a short 
distance between exposures, then the diffraction halo generated by the 
processed plate will consist of a pattern of parallel straight fringes 
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similar to Young's fringes. The diffraction halo observed through 
a small area of the recorded image will correspond to the local dis- 
placement at the corresponding point on the object and the direction 
of the fringes will be orthogonal to the direction of the local dis- 
placement vector. Also, Burch and Tokarski [5] showed that by optically 
illuminating the developed photographic plate with a converging spheri- 
cal wave, the entire surface could be analyzed at one time to determine 
the displacement field of the surface. 

Leendertz [6] developed a technique utilizing laser speckle effect 
for measuring either normal or in-plane components of displacement over 
an entire surface at one time. For measurement of normal (out-of-plane) 
displacement two surfaces are coherently illuminated and two separate 
speckle patterns are obtained (Figure 1). For measurement of the in- 
plane components of displacement, a surface is illuminated by two beams 
of coherent laser light, symmetrically disposed about the normal to the 
surface as illustrated schematically in Figure 2. These two speckle 
patterns are superimposed and the resultant speckle pattern is recorded 
on film. The intensity distribution of the resultant speckle pattern 
depends on the relative phases of the component patterns. Then one or 
both speckle patterns is changed and again the resultant speckle pattern 
is recorded on the same photographic film. By measuring correlation be- 
tween the resultant pattern at two different times, a change of relative 
phase is detected, which in turn gives a measure of surface displacements. 
These correlation fringes are observed either in real-time or by combining 
two transparences having resultant speckle pattern at two different times 
and illuminating it in a Fourier filter system. A major drawback of 
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this technique is the path length difference between the two illumi- 
nating beams has to be less than the coherence length of the light 
used to generate correlation fringes. 

Duffy [7] presented two methods for measuring in-plane surface 
displacements. In the first method, this displacement was determined 
by photographing a coherently illuminated object through two laterally 
displaced apertures. The displacement is displayed as a pattern of 
Moire 1 fringes over the image of the surface. Thus there is no need 
for scanning of the beam on a point by point basis. As the surface is 
illuminated by only a single laser beam, the implementation problems 
associated with the dual-beam technique (mechanical stability and 
equal path lengths between the various optical components) are mini- 
mized. In the second method, the object is illuminated using a single 
laser beam and photographed via a double exposure before and after dis- 
placement. The Fourier transform of the doubly exposed transparency 
was obtained optically by illuminating the photographic plate with a 
converging spherical wave. The main advantage of this procedure is that 
the whole-field displacement can be analyzed and by appropriate position 
of a set of apertures in the transform plane any component of the dis- 
placement normal to the line of sight can be detected and with variable 
sensitivity. This thesis presents the mathematical analysis which 
describes the formation of fringes for both the whole field and point- 
by-point method of data reduction. 

In addition, research has recently been conducted into the use of 
laser speckle interferometry [8, 9, 10, 11] for vibrational analysis. 

In this thesis, the time-average theory using the Fourier transform 
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is developed to present the application of this technique to measure- 
ment of in-plane displacement induced by the vibration of an object. 




'll. DOUBLE EXPOSURE SPECKLE INTERFEROMETRY 
2 . 1 Basic Phenom ena 

When using a coherent laser light Lor illumination, a uniformly 
diffuse surface has a speckle or grainy appearance due to random in 
terference within the resolution of the eye (or photographic system). 
Additionally, any point in front of this surface will receive contri- 
butions from all points on the surface which have similar amplitude 
but random phases. Thus the amplitude and phase of the radiation 
field in front of the object vary in a random manner from point to 
point with the average intensity being the same as if the light were 
incoherent. This phenomenon is known as the speckle effect and has 
been investigated by several authors [12]. When an- imaging system 
is used, the size of the speckle in the image plane is inversely 
proportional to the aperture because intensity variations cannot be 
produced inside distances less than the Airy disk diameter [13]. The 
speckle appearance of an illuminated object and the effect of aperture 
on speckle size is illustrated in Figure 3. 

The speckle effect provides a sensitive method for measuring dis- 
placements and strains on the surface of an object. The correlation 
between the relative phases of speckle patterns of a coherently illu- 
minated object, at two different times, is the basic principle of the 
Leendertz technique [6]. With this technique, a double exposure 
specklegram (photo-negative) is made of the illuminated surface before 
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and after deformation. Thus, each bright "speckle point" is recorded 
on the specklegram as two dark points: one in the undisplaced position 

and the other in the displaced position. Finally, only knowledge of 
the photographic magnification factor is required to determine the 
local in-plane displacement once the length and direction of the line 
segment between two speckle points are known. 

In this report, the major emphasis is placed upon the "whole-field" 
interpretation of Young's fringes and comparison of the results to the 
"point by point" interpretation. The theory presented here extends 
works conducted by Burch and Tokarski [5] and Duffy [7] in the inter- 
pretation and data collection of the fringes. 

One method for determing displacements from the specklegram, as 
presented by Kinariwala- [14], is to direct a narrow collimated laser 
beam through the specklegram. Diffraction will modify the emerging 
light rays into a cone. This "diffraction halo" is the result of 
diffraction from the random distribution of small speckles. Since 
the speckles are recorded in pairs (displaced and undisplaced), 
a parallel fringe pattern (Young's fringes) also occur in the 
"diffraction halo", (Figure 4). From these fringes, displacement 
at a point can be determined. Thus, the entire in-plane displacement 
field can be determined by mapping the surface on a point-by-point 
basis. The number of points analyzed is dependent upon the type of 
analysis conducted. To decrease the uncertainty in the displacement 
field, more points would have to be analyzed. 

In the second method of data interpretation, a converging spherical 
wave of coherent light is passed through the specklegram resulting in 


Direction or inpiane 
Displacement 



Figure 4. Arrangement for Analysis of the Speck! egram - Point-by-Point Interpretation 
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the specklegram being imaged in the transform plane as shown in Figure 
5. Using this "whole field" technique, the displacement field of the 
entire surface in any direction can' be determined with a single 
photograph. The sensitivity of the displacement magnitude is 
dependent upon the amount of offset of the aperture, p, in the 
transformed plane; while the displacement component orientation 
depends upon the angle between the focal point in the transform plane 
and the aperture. For example, with a zero angle between the two, 
only horizontal in-plane motion would be measured. While with a 
90° angle, only vertical in-plane motion would be measured. 

2.2 Theory of Speckle Interferometry - Double Exposure - Single 

Beam Analysis 

Consider the arrangement shown in Figure 6. An object is 
illuminated by a single laser beam {monochromatic and coherent). 

The location of the light source is shown as S which illuminates 
a pon'int P on the surface of the abject. In this analysis, the 
"double-exposure" technique for obtaining the specklegram is used, 
meaning that initially an exposure of the undeformed surface is taken,' 
then the surface is deformed and another exposure is taken on the 
photographic plate. In this development, the intensity distribution 
on the photographic plate is investigated and, using the results of 
this analysis, the in-plane displacement of the exposed surface is 
determined. 

From the first exposure, the complex light amplitude can be ex- 
pressed as 

E (x-j sXg) -A (x-j jXg) exp £i® (x-j,x^)j 


( 2 . 1 ) 
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where 

(x-j, x 2 ) - Film plane coordinates 

0 (x-j, x 2 ) - Relative phase at each location cn the film 
A (x-j, x 2 ) - Amplitude factor 
The intensity of the first exposure is given by 

Ij = Ej E* (2.2) 

where is the complex conjugate of E-j 

I-j = A (x-j , x 2 ) exp L^(x 1 , x 2 )] A (x-j , x 2 ) exp C-iefx^ x 2 )] 

(2.3) 

1 1 = A 2 (x r x 2 ) (2.4) 

For the second exporure, each point on the body is displaced to 
a new position; therefore, the exposure is 

E 2 (x^, x£) = A (x^, x£) exp [i0(xj, x£)] (2.5) 

and the intensity is 

1 2 = A 2 (x-j , x 2 ) (2.6) 

Thus, the total intensity on the photographic plate is the total of I-j 
and In. 

C I-J- “ I-J + I 2 (2.7) 

ly = A 2 (x-j , x 2 ) + A 2 (x-j , x£) (2.8) 

What is needed is to relate this intensity distribution to the 
displacement of the surface. In the first exposure, the point P is at 
coordinate (x-j, x 2 ). For the second exposure, the point P is displaced 
to point Q and is imaged on the film plane at coordinate (x-j, x 2 ). 

The displacement vector on the surface of the object, represented by 
PQ, is 

pq = u,?, + u 2 i 2 


(2.9) 
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and the displacement vector in the film plane is 

^Film Plane = u lf e l * u 2f e 2 ^ 2 * 

and the displacement in the film plane in relation to the object is 


If = Hu l 


(2.11-1) 

2f ‘ ^ u 2 

■» 

(2.11-2) 


where M is the film magnification factor, usually less than 1. 

Therefore, it follows that the film plane coordinates for the second 
exposure is given by 

x-j = x-| + u-| f * (2.12-1) 

x 2 = x 2 + u 2f (2.12-2) 

With these two relations, the total intensity given in Equation (2.8) 
can be rewritten as 

I-j- = A (x-j , X 2 ) t A (x-j + u-j x 2 t Ugf) (2.13) 


If the exposure E(x-j, x 2 ) at each point of the plate represents 
small fluctuations which are mostly confined to the toe and the bottom 
portion of the H and D curve, the amplitude transmission g(x-j, x 2 ) of 
the photographic plate is approximately a linear function of the exposure 
i.e. 


g(x ls x z ) = a + b I T 


(2.14) 


where a and b are constants of the photographic film. With the result 
in eq (2.13) this can be written as 

gfxp x 2 ) = a + b [A ( X-j , 'X 2 ) t A (x-j ■** u-j^, x 2 + u 2 ^)j 


(2.15) 
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This expression gives the relationship between transmissivity and 
intensity in terms of the displacements in the film plane. The next 
step in the analysis is to relate the displacement in the film plane 
to some measurable quantities. In this development, the "whole field" 
interpretation of the fringes was used in the determination of dis- 
pl acement. 

2.3 Optical Fourier Fringe Interpretation 

- The whole field fringes are obtained by taking optically the 
Pourier transform of the amplitude transmission function g(x^, Xg) of 
the specklegram [5]. To simplify the notation, the Xg coordinate will 
be suppressed in the following development, i.e. 

g(x) = a + b [A 2 (x) + A 2 (x + u f )] (2.16) 

i t 

'This optical transform is accomplished by illuminating the specklegram 
with a converging spherical wave as shown in Figure 5. Neglecting 
aberrations, the light amplitude G(w) in the transform plane is propor- 
tional to the Fourier transform of g(x) times a quadratic phase factor 
[5] as given by the following equation 

t 

. 2 

G(w) = exp (-g-jp- ) ./’gC x) exp (-iwx) dx (2.17-1) 

where 

w - k£_ - k tan 0 (2.17-2) 

2 

where p is the coordinate of the iris in the transform plane and 

k = 2lT / X 5 where x is the wave length of the converging beam and 

0 is defined in Figure 5. The direction of the laser beam and the 
common axis of the lenses L-j and are assumed perpendicular to the 
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transform plane and to the plane of the specklegram. Equation (2.17-1) 

follows from the Kirchoff diffraction formula which was derived by 

Burch and Tokarski [5], Eq (2.16) can be rewritten as 

. 2 

G(w) = exp [ |p ] /{a + b [A (x) + A (x + u f )]} exp (-iwx)dx 

(2.18) 

.2 .2 

G(w) = a exp [ |p ] /exp (-iwx)dx + b exp [ |P ] 

S [A 2 (x) + A 2 (x + u f )] exp (-iwx)dx (2.19) 

The first integral in eq (2.19) is the delta function, <s(w); therefore, 

6(w) = a exp [ ] s(w) + b exp [ |3T ] / [A 2 ( x ) + 

A 2 (x + u^)] exp (-iwx)dx (2.20) 

Equation (2.20) can be written in the following form with the applica- 
tion of the shift theorem. 

- 2 .2 

G(w) = a exp [ |p ] s(w) + b exp [ |p ] D + exp (-iwu f )] 

f A 2 (x) exp (-iwx)dx (2.21) 

The last integral is the definition of the Fourier Transform of A (x) 

2 

and is represented as F [A ] which leads to 

2 2 

- G(w) = a exp [ Ift ] <(w) + b exp [ if- ] F [A 2 ] 

[1 + exp (-iwu^)] e (2.22) 

The delta function represents the (idealized) point focus of the 
illuminating beam. It contributes to the diffracted amplitude only 
at the point w=0; in practice, it is a small area around the point. 
Outside this small area, the delta function is zero and G(w) is given by 

G(w) = b exp [ |p ] F [A 2 ] [1 + exp (-iwu f )] 


(2.23) 
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The next step in the development is to determine the intensity 


function in the transformed plane which is 
X F = G(w) G*(w) 

where G*(w) is the complex conjugate of G(w) 

I F = ( b exp [ ] F [A 2 ] [1 + exp (-1wu f )3 ) 

C b exp [ ^jr 3 F [A 2 ] ci + exp (1wu f )] ) 

Ip ■ b 2 F [A 2 ] 2 [2 + exp (-1wu f )+ exp ( l'wu,) ] 
Using the tri gome trie identity 


(2.24) 


(2.25) 


(2.26) 


cos (wup) = exp (iwUp) + exp (-iwup) 

2 

the final expression for the intensity is 
I F = 2b 2 F[A 2 ] (1 + cos wu f ) 


(2.27) 


(2.28) 


A fringe will be defined when the intensity. Ip, is equal to zero. 
2 

Since b and F [A ] are both nonzero, fringes will occur when 


[1 + COS WUx] = 0 


(2.29) 


' COS WUp = -1 


(2.30) 


which means 


wup = (2n - l)ir n = 1, 2, 3, 


(2.31) 


With the use of eq (2.17-2), eq (2.31) can be put in the following form 


2-frpu.c 

— r— L = (2n - 1) it 


(2.32) 


U P " (n - j) p 


(2.33) 


Eq (2.33) expresses the in-plane displacement in the film plane in terms 
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of quantities which can be measured in the transform plane. Finally, 
the displacement of the deformed surface is given by 

u « ^f (2.34) 

m 

“ = <"-?>< sf> < 2 - 35 > 

Information at location p is selected by placing an aperture to 
pass only the light in the immediate neighborhood of p as shown in 
Figure 7. Cloud [12] reported that an aperture opening of approximately 
one-tenth the amount by which the aperture is displaced from the optical 
axis gives the best results. This light is then imaged on film. Thus 
p being fixed, the image fringe information is proportional to u^ as 
shown in eq (2.33). The direction of the measured displacement is 
determined by the relation of p with the u- and v- axis as shown in 
Figure 7. In order to measure the horizontal displacement, the 
aperture is moved along the u-axis and to determine the vertical dis- 
placement, the iris is moved along the v-axis with movement always 
starting from the focal point in the transformed plane. Thus, in- 
plane displacement in any direction can be determined by appropriate 
selection of the aperture coordinates. Eq (2.33) illustrates that 
the sensitivity of the displacement determination is based on the 
distance of the aperture from the focal point, p. An increase in 
p causes an increase in the sensitivity of the measurement; however, 
due to the resolution limit of the film and the relative width of the 
fringes, the sensitivity of this technique is limited, dependent upon 
the amount of total displacement. 


CmJ 



Figure 7. Arrangement for Fringe Interpretation - Whole Field Technique 
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2.4 Point-by- Point Data Interpretation 

The point-by-point data interpretation can also be determined 
from eq (2.33). Figure 8 is a schematic of the point-by-point data 
analyzer. A narrow beam of coherent light is passed through a region 
of the film and fringes occur in the observation plane. To correctly 
interprete these fringes, a reexamination of eq (2.33) is needed. In 
the development of the transformed light amplitude (eq 2.17-1), inte- 
gration over the entire film plane was conducted. For the point-by- 
point analysis, this intensity relation is zero except in the narrow 
region of the film through which the light passes. The in-plane film 
displacement is given by eq (2.33) 

u f = (n - \ > r 

~ For this case, the in-plane displacement is constant with p being the 
variable. From the worked conducted by Kinariwala [14], the expres- 
sion for the in-plane film displacement for n = 1 



where 

A - wavelength length of light source 
z - distance from film plane to observation plane 
S - distance between fringes 
For the case of n = 1, eq (2.33) 


u 


f 


„ Jz 
“Zp 


(2.36) 


(2.37) 


-An examination of the last two equations, shows the displacement 
determined from eq (2.36) to be twice eq (2.37). However, this 
difference arises in the measurement of X and p. The distance X 
is measured between the fringes while p is measure from the center of 


the halo to the nearest fringe. Figure 8 presents these two measure- 
ments. Thus, for the case of n = 1, 

X = 2p (2.38) 

and the two equations are equivalent. 

2.5 Experimental Verifications 

Results of experimental verifications of the development in 
section 2.3 are reported in this section. Two examples are considered 
and for each case theoretical and point-by-point results are presented 
for comparison to the whole field results. An account of experimental 
procedures of recording and collection of data is also presented. 

Experimental Equipment 

The experimental set up for recording a specklegram is shown in 
Figure 9 and schematically in Figure 10. An air-supported table was 
used to isolate the optical elements from surroundings and to provide 
a vibration free support system. 

A Spectra- Physics model 125A He-Ne cw laser operating at 80 mw 
plane polarized output was used as the light source. The light was 

spatially filtered using a Spectra- Physics spatial filter 332. Light 

reflected from the object was imaged onto Kodak High Speed Holographic 

Film 131-02 using a Wollensok 15 inch f/5.6 lens. The emulsion was 

supported on a 4 in. x 5 in. glass plate to prevent reduction of the 

imaged picture. 

In the point-by-point fringe interpretation, a set up similar to 
the one shown in Figure 8 was used. The light source was a Spectra- 
Physics Stabilite™ Model 120 He-Ne cw laser. The diffraction halo 
with fringes was photographed with a 35mm Nikon camera using Kodak 
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C - Camera 


Figure 9. Experimental Set up for Taking a Specklegram 
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Figure 11 
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Panatomic-X film. For the whole-field fringe interpretation, the 
set up shown in Figure 11 was used. The Spectra-Physics Model 125A 
was the light source. The two lenses were spherical lenses with a 
focal length of 26.5 inches. The fringes were photographed using 
a 35mm Nikon camera with Kodak Panatomic-X film. 

In the first example, only the vertical displacements were de- 
termined for comparison to the theoretical results. For the second 
example, horizontal displacements were determined. 

Example 1 - Cantilever Beam 

The first experimental example chosen to verify the theoretical 
results was the bending of a rectangular cantilever beam with an end 
load, as shown in Figure 12. The beam had the end fixed at ^ 0 
and the end z 1 = 1 had a single force P directed along the positive x- 
axis producing a known deflection, 6 , at z 1 = 1. 

The z^-axis was taken along the center line of the beam and the 
x- and y-axes were orthogonal axes intersecting at the centroid of 
the fixed end. The vertical (u) displacement component as given by 
Sokol nikoff [15] is listed below 



The load P is governed, by the deflection of the free end. Suppose 
the load P is applied at ^ 1 such that u(0,0,l) * 6 where 6 is the 
known deflection. From eq (2.39) 

P = (2.40) 

r 

U = (1-z^ (x 2 - /) - f 3 - If'] 


*2.41) 



29 


Material used for the beam was plexiglass and the physical pro- 
perties of the beam are 
1 =6.00 inches 
b = 0.125 inches 
c = .255 inches 
v = .37 
E = 456000 psi 
Si = .003 inches 
$2 = -006 inches 

Fixed end conditions for the beam were obtained by bonding one 
end of the plexiglass to an aluminum block as shown in Figure 13. The 
end load was applied with a micrometer screw and the deflection was re- 
corded with a dial indicator (Figure 14). Deflections of .003 inches 
and .006 inches w^re investigated. 

Discussion and Results of Cantilever Beam Experiment 

Figure 15 presents the experimental set up. The coordinates rela- 
tive to the cantilever beam are given in Figure 12. Theoretical results 
for 6 = .003 inches and 6 = .006 inches are tabulated in Table 3^ 
and Table 4, respectively. Also, a plot of deflection versus distance 
along the beam for the two loading conditions are shown in Figure 16 
and 17. 

As discussed in section 2.4, using the point-by-point data reduction 
method, the in-plane displacement in the principal plane is given by 


u 


- hl. 
f ~ 2P 


(2.42) 


Figure 18 shows the diffraction halo and fringes at various locations 


Figure 13. Fixed end Condition for Cantilever Beam 



Figure 14, Method of Load Application and Measurement of <5 







32 


along the beam for <5 = .003 inches. The in-plane vertical displace- 
ment is given by 

u = u_ sin 6 {2.43} 

f 


wnere 6 is measured as shown in Figure 3. The data taken for the two 
loading conditions are tabulated in Table 1. For comparison purposes, 
the experimentally determined displacements are plotted on Figure 16 
for 5 = .003 in. and Figure 17 for 6 = .006 in. 

For the whole-field fringe interpretation, the in-plane displace- 
ment in the film plane was shown to be equal to 

u f ■ (n - r> t (2 - 33) 


and the actual displacement is 

- ( n h Az 
j (n - qJ 

2' mp 


(2.35) 


Figure 19 is a photograph of the fringes in the transformed plane for 
<5 = .006 inches. For these two photographs, 

p = .98 inches 
z = 33 inches 
x = 20.22 microinches 

The data for both load conditions is tabulated in Table 2. Again, 
for comparison, the experimental results are plotted on Figure 16 for 
6 = .003 in. and on Figure 17 for 6 = .006 in. 

As seen from the results presented in Figures 16 and 17, excellent 


results were obtained for both techniques. The maximum difference in the 
point-by-point interpretation was ^5%. For the whole-field interpre- 
tation, the maximum error was also ^5%. 
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TABLE 1 

EXAMPLE 1 - CANTILEVER BEAM - POINT-BY-POINT INTERPRETATION 
FOR VERTICAL DISPLACEMENT 


- .003 inches 


= .006 inches 


z 

P 

0 

u 

in) 

'(in) 

•Odea) 

(in) 

0 

.148 

- 3.0 

2.88 x 10' 

. .3 

.148 

- 3.0 

2.88 

.6 

.148 

- 3.0 

2.88 

.9 

.177 

- 3.0 

2.40 

1.2 

.197 

- 4.0 

2.15 

1.5 

.216 

- 4.0 

1.95 

1.8 

.241 

- 4.0 

1.75 

2.1 

.276 

- 4.0 

1 .54 

2.4 

.315 

- 6.5 

1.34 

2.7 

.364 

- 6.5 

1.16 

3.0 

.433 

- 6.5 

.98 

3.3 

.521 

-10.5 

.80 

3.6 

.630 

-14.0 

.65 

3.9 

.748 

-15.0 

.51 

4.2 

1.164 

-18.0 

.39 


Z 

ini 

P 

(ini 

6 

f dec) 

u 

. ( ini 

0 - - 

.3 

. .375 

- 3.0 

5.27 x 10 

.6 

.400 

- 3.0 

4.94 

.9 

.450 

- 3.5 

4.40 

1.2 

.475 

- 3.5 

4.16 

1.5 

.500 

- 3.5 

3.96 

1.8 

.625 

- 4.5 

3.44 

2.1 

.675 

- 4.5 

2.93 

2.4 

.800 

- 7.0 

2.48 

2.7 

.925 

- 7.0 

2.14 

3.0 

1.025 

- 7.0 

1.92 

3.3 

1.200 

- 7.0 

1.64 

3.6 

1.425 

-7.5 

7.28 

3.9 

1.765 

-12.5 

1 .00 

4.2 

2.325 

-12.5 

.76 








TABLE 2 


EXAMPLE 1 


- CANTILEVER BEAM - WHOLE FIELD INTER 


FOR 

VERTICAL 

DISPLACEMENT 


<5 = .003 

inches 

n 

2 

u 

— 

(inch) 

(inch) 

1 

3.86 

.487 x 10" 3 

2 

2.13 

1.461 

3 

.77 

2.436 


6 = .006 

inches 

n 

z 

u 

— 

(inch) 

( inch) 

1 

4.63 

.487 x 10' 3 

2 

3.36 

1.461 

3 

2.45 

2.436 

4 

1.72 

3.410 

5 

1.09 

4.386 

6 

.36 

5.360 
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TABLE 3 

. ... EXAMPLE 1 - CANTILEVER BEAM - THEORETICAL RESULTS 

FOR VERTICAL DISPLACEMENT 







Vertical Displacement (inch x 10) 


* 
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Example 2 - Pressure Vessel 

The second experimental example chosen to verify the results 
of Section 2.3 was an internally loaded thin walled cylindrical 
pressure vessel. The coordinate system for this example is shown in 
Figure 21. The y-axis was chosen at the axis of symmetry. With 
the x- and z-axes intersecting at the center of the base of the 
cylinder. The z-axes on the model and on the film were parallel. 

The radial strain for a capped end thin wall pressure vessel 


is 


where 


A R 1 / \ 

£r = R “ ~ E °0 ~ 


„ _APR 
0“ t 


„ - APR 

z “ 2t 

Therefore, the radial displacement is 

AR-lf 2 <2-v) 

The horizontal displacement, u, is 
u = A r sin Q 

However, 

sin 0 =4 

Thus, the theoretical horizontal displacement is 
u = (2-v) 


(2.42) 

(2.43-1) 

(2.43-2) 

(2.44) 

(2.45) 

(2.46) 

(2.47) 


For this example 
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AP « 400 psi 
R = 2.54 inches 
t = .09 inches 
E - 10 x 10 6 psi 
v = .33 

A tabulation of theoretical horizontal displacement for various x's 
is presented in Table 6. 

The capped end conditions for the pressure vessel are obtained by 
welding aluminum plugs at both ends. The physical dimensions for the 
pressure vessel are shown in Figure 22. 

The same procedure was followed as in the first example. The 
results for the point-by-point interpretation are tabulated for example 
2 in Table 4. Results for the whole field interpretation presented 
in Table 5. A photograph of the whole field fringes is shown in Figure 
23. 

To verify the speckle theory, a comparison of the experimental 
results to the theoretical results is presented in Figure 24. Excellent 
agreement is obtained near the center of the cylinder; however, once 
removed from the center area, a small difference between the theory and 
experimental results occurs. For the whole field interpretation, rela- 
tionship between x and the horizontal displacement is a straight line 
as would be expected; however, the slope of the curve is less than 
the theoretical. For the point-by-point results, initially the rela- 
tionship between x and the displacement component is linear; however, 
at larger values of x, this linearity disappears. Probable cause for 
this difference is in the focusing of the cylindrical body onto the film 
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plate. If the center of the cylinder is used for the point of focus, 
the outer edges (as x increases) will be slightly out of focus depen- 
dent upon the radius of the cylinder. Recent work has been conducted 
into this problem of misfocusing and has shown that this can lead to 
a difference in fringe spacing. 


TABLE 4 


EXAMPLE 2 - PRESSURE VESSEL - POINT-BY-POINT INTERPRETATION 
FOR HORIZONTAL DISPLACEMENT 


X 

P 

0 

u 

(inch) 

(inch) 

(deq) 

(inch) 

0 

.25 

0.0 


.2 

.24 

7.5 

1.93 x 10' 4 

.4 

.23 

15.0 

3.91 

,6 

.22 

21 .0 

5.92 

.8 

.21 

28.0 

7.94 

1.0 

.20 

31.0 

8.92 

1.2 

.19 

36.0 

10.44 


TABLE 5 

EXAMPLE 2 - PRESSURE VESSEL - WHOLE-FIELD INTERPRETATION 
FOR HORIZONTAL DISPLACEMENT 


X 

(inch) 

n 

u 

(inch) 

.17 

i 

1.619 x !0‘ 4 

.52 

2 

4.858 

.88 

3 

8.096 

1.25 

4 

11.330 

1 .66 

5 

14.570 




TABLE 6 


EXAMPLE 2 - PRESSURE VESSEL - THEORETICAL RESULTS 
FOR HORIZONTAL DISPLACEMENT 


X 

u 

(inch) 

(inch) 

0 

0.00 

.2 

1.88 x 10“ 4 

.4 

3.77 

.6 

5.65 

.8 

7.54 

1.0 

9.43 

1.2 

11 .31 

1.4 

13.21 

1.6 

15.14 

1.8 

17.02 
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I 

| Fiqure 24. Frinqes from Whole Field Interpretation for Pressure 

Vessel, aP = 400 psi 

1 

I 


III. TIME AVERAGE SPECKLE INTERFEROMETRY 


3.1 Basic Phenomena 

As presented in Section II, in-plane displacement resulting from 
a static load can be determined using double exposure speckle inter- 
ferometry. In this section, the speckle theory is extended to in- 
clude the time averaging technique for use in determining in-plane 
displacement due to harmonic vibration of an object. The time averag-, 
ing technique has previously been sucessfully used in holography [16]. 

In recent years, time averaging in conjunction with speckle interfero- 
metry has been investigated. Tiziani [10] showed the intensity expres- 
sion in the transformed plane to be a function of the zero order Bessel 
function for a time average recording of a harmonically oscillating 
object or image. The work conducted in this section of the report 
extends the work conducted by Tiziani [10] by determining the relation- 
ship between the transformed intensity function and the in-plane 
displacement of the oscillating object. 

The analytical work closely parallels the development of the double 
exposure technique presented in Section II. The major difference is 
in the determination of the intensity expression. For the double 
exposure technique, the intensity is determined for each exposure and 
summed. For the time averaging technique, integration is performed 
over the total time to arrive at the intensity expression. 
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3.2 Theory of Speckle Interferometry - Time Averaging - 
Single Beam Analysis 

The arrangement presented in Figure 6 is also used in the time 
averaging. An object is illuminated by a single laser beam with the 
location of the light source shown as S which illuminates a point P 
on the surface of the object. In time averaging, a single exposure 
is taken of the vibrating surface with an exposure time at least 
equal to the time for one complete cycle of the vibration. 

For the condition of the time dependent loading, the complex 
light amplitude can be expressed as 

E (x jj x 2 ) — A (Xj t x p x 2 t x^) 

exp i[ 0 (x ls x 2 ) +A0] (3.1) 

where 

x i’ x 2 " P^ ane coordinates 

x j, Xg - Function of the amplitude of vibration, 
frequency, time and location 

0 - Phase angle of the beam 

A0 - The relative change in the phase angle 
due to the vibration 

For the time averaging, the film exposure over a time interval 

t is 

tl y = J Etx i# x 2 ) E* (Xj, x 2 ) dt (3.2) 

O 

where F* (x lS x 2 ) is the complex conjugate of F (x , x 2 ). Again, to 
simplify the notation the x 2 -coordinate will be suppressed for con- 
venience, i.e. 

A _ 

tI T = E(x) E*(x) 


dt 


(3.3) 
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Therefore . 


tl.j. = A (> 


x + x 1 ) dt 


93.4) 


A physical interpretation of x 1 is needed before this development 
proceeds. At t = 0, the point P is at coordinate (Xj, X 2 ) on the object. 
For t > 0, the point is displaced to a new location P which is dependent 
upon the amplitude of vibration, the frequency and the time. Thus, the 


line segment PQ can be written as 


PQ = u (X x , X 2 ) sin tut (3.5) 

where u represents the maximum displacement of the surface at (X x , X 2 ) 


and on the film plane 


™ Film Plane = u f {x i’ sin * 


(3.6) 


u f 

u = r 


where M is a magnification factor. Thus 


(3.7) 


tl t = j $ (x + sin tut) 


(3.8) 


As present in Section 2.3, the amplitude transmission g(x) can 


be approximated by a linear function of intensity, i.e, 
g(x) = a + b I T 

where a and b are constants of the film. Thus, 
g(x) = a + 6- / t A 2 (x + Uf sin to t) dt 


(3.9) 


(3.10) 


Eq (3.10) is a relationship between transmissivity and intensity in 
terms of the displacement in the film plane. 
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3.3 Optical Fourier Fringe Interpretation 

The whole field fringes are obtained by taking optically the 
Fourier transform of the amplitude transmission function g(x^» Xg) 
as discussed in Section 2.3. The light amplitude function G(w) *n 
the transformed plane is [5] 

G(w) = exp ) f g(x) exp (-iwx) dx (3.11) 

as presented in Section 2.3. 

For the time averaging, 

2 

G(w) = exp (|j~ JJ^a + £ [ A 2 (x + u f sin tot) dt 

exp (-iwx) dx (3.12) 


or 

G(w) = exp (^j— ) J* a exp (-iwx) dx + b (lip. ) 


//^ 2 ( x + u f sin wt ) ex P (~i wx ) dt dx 

o T 

(3.13) 

Now 


6(w) = f exp (-iwx) dx 

(3.14) 

Thus 

G(w) = afi(w) exp (^jp ) + | exp (^jp ) 



(x + u f sin tot) exp (-iwx) dt dx (3.15) 

0 T 

2 

Since t and x are independent, and A (x + u^ sin tot) and exp (-iwx) 

are continuous functions-, and u^, is a function of x only, the second 

part of eq (3.15) can be written as 
2 

T exp ^ 'C'f A 2 (x + u^sin tot) exp (-iwx) dx dt 


(3.16) 
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From eq (2.21), the first integration is the Fourier transform of 
2 

A (x + u.pSin wt). Thus, 

G(w) = as(w) exp ) + ^ ex P ) 

F [A 2 (x.j + u^r sin tot)] dt (3.17) 

Employing the shift theorem, eq (3.17) can be rewritten as 
G(w) . as(w) exp (|f 2 ) ♦ exp (i™ 2 ) 

f^exp (iUxW^in tot) dt (3.18) 

Jo 

As discussed in Section 2.3, the delta function's contribution occurs 

only at w = 0„ For any value of w other than zero, eq (3.18) is 

2 2 

G(w) = bF ^ (x - ^ - exp )£ exp (iu f wsintut) dt (3.19) 


From Bessel function analysis 

2n0 o {x) =^ 2ir exp (ixcosQ) dQ =£ 2ir exp (ix sine ) de (3.20) 
where J Q (x) is the zero order Bessel function. For evaluation of 
eq (3.19), the integration will be for one cycle, that is t will be 
the time required for one cycle, thus if many cycles occur the time 
will be some multiple of t. Rearranging eq (3.19), and multiplying 
and dividing by to gives 

G(w) = tu exp ^ F ^M^/J’exp ( i u f w sin “t) tudt 

(3.21) 


Comparing eq (3.21) with eq (3.20) 


q = ait 


de = oidt 
x = u f w 
t = 2 tt/w 
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Consequently, eq (3.21) becomes 

S(w) - ^ exp ] F [A 2 ] J. (u f w) (3.22) 

The intensity function in the transformed plane is 

Ip = G(w) G*{w) (3,23) 

Therefore, 

Ip = -fy F [A 2 ] 2 j! (u f w) (3.24) 

t ^ 

Fringes are defined in the transformed plane when Ip = 0. Since all 
the terms except J 0 (u f w) are always non-zero, fringes occur when 

Jo (u f w) = 0 (3.25) 

From eq (2.17) 


w ® 


2irp 

XZ 


Thus 


■2ttPu.p 

o.(^) - 0 


(3.26) 


Since all the terms inside the bracket except u^ can be measured or 
are known, the in-plane displacement for a harmonically oscillating , 
object can be determined. Table 7 is a tabulation of the zero 
values for J 0 (x) . 

' 3.4 Point-by-Point Interpretation 

Similarly to the discussion in Section 2.4, the point-by-point 
interpretation can be determined from eq (3.26). For the point-by- 
point interpretation, the in-plane displacement is constant at each 
location with the distance p being variable. For the case of n - 1, 
eq (3.26) yields 
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TABLE 7 

ZERO VALUES FOR J 0 (x) 
J 0 (x) Equals Zero when 


X 01 

x 02 

X 03 

x 04 

x 05 

X 06 

X 07 

x 08 

x 09 

X 10 

X 11 

x 12 

X 13 

x 14 


= 2.40 
= 5.52 
= 8.65 
= 11.79 
= 14.93 
= 18.07 
= 21.21 
= 24.34 
= 27.49 
= 30.63 
= 33.78 
= 36.92 
= 40.06 
= 43.20 
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2irpu i 

~Tz~ 


2.40 


(3.27) 


or 


.. _ .382XZ 

u f ' p 


(3.28) 


and 


u f = 


.382Xz 

Mp 


(3.29) 


Comparison of eq (3.28) to eq (2.37), shows the only difference 
between the time averaging and double exposure point-by-point inter- 
pretation is a difference in the constant. 

3.5 Numerical Example - Vibrating Cantilever Beam 
The numerical example chosen is a vibrating cantilever beam. 

The beam is made of aluminum. and the dimensions of the beam are pre- 
sented in Figure 25. For this example, the first two modes of vibra- 
tion were analyzed. 

The solution for the amplitude W(z^), as presented by Nowacki 
[18] for free transverse vibration of a cantilever beam with the end 
fixed at z = 0 and an initial tip deflection of 5 at z = 1 is for 
each mode of vibration 

S( a. ) 

(3.30) 


s(e ) 

W r (2 1 ) - C CU(x r z 1 ) - r ° 1<z 


where C is dependent upon the initial conditions and 


utvj) - 


cosh LZj - COS X y ,Z 1 


VU^) = sinh X r z^ - sin X r z^ 


2 
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cosh 8 + COS 8 

s <«r> r -2 E 

sinh 8 + sin 

T(B r ) 4 11 



Bj = 1 .875 

3 = 4.694 
2 

- natural frequency 

If the initial condition is such that only the first mode is 
activated, then at z =1 

H t (i) -« - c [utBj) - tt^t v ( e j.)] (3 - 31) 

Therefore, for this case 


6 


S(B ) 

c [“(A 2 !* - nfr vtJ i z i )] 
sTsTT 

c cu(s i> - TTe^T v(8 1 )] 


(3.32) 


Similarly, if the initial condition is such that only the second mode 
is activated then 


S(BJ 

W 2 ( Zl ) [U( ^ 2 z i ) - v( \ z i)] 

~ = sTbT] 

CU{3 2 ) -T^yV^)] 


Therefore, for each case, the ratio of the amplitude to the initial 
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tip deflection as a function of distance along the beam, z 1# can be 
determined and is tabulated in Table 8. 

The amplitude of vibration is just the in-plane displacement, 
u, defined in Section 3.2, i.e. 

For the first mode of vibration 


S(Bi) 

[U{\ 1 z 1 ) y V(? 1 z 1 )] 

5(0 

" T^r V(S i )] 


(3.34) 


The in-plane displacement in the film plane is 

Uy = Mu (3.35) 

where M is a magnification factor. 

Then for the first mode of vibration, the in-plane displacement in the 
film plane is 

s(e x ) 

[UfAjZj) ^ v( Vl )] 

u f = M SC] < 3 - 36 > 

" T(0 X ) 

For the second mode 


S(3 2 ) 

[U(*2 Z l) “ T(0 ) V(A 2 Zi)3 

U f - M — (3 * 37j 

[u(p 2 ) - y^j v(e 2 )] 

Whole Field Interpretation 
For the whole field interpretation, let 


p = .3 inches 
z = 26.5 inches 
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TABLE 8 

RATIO OF AMPLITUDE/TIP DEFLECTION 
FIRST MODE SECOND MODE 


Z 1 

u/s 

Z 1 

ll/S 

inch ) 

— — 

( inch ) 


0 

0 

0 

0 

.3 

.0042 

.3 

-.0254 

.6 

.0165 

.6 

-.0926 

.9 

.0363 

.9 

-.1888 

1.2 

.0630 

1.2 

-.3011 

1.5 

.0961 

1.5 

-.4173 

1.8 

.1349 

1.8 

-.5262 

2.1 

.1789 

2.1 

-.6178 

2.4 

.2275 

2.4 

-.6836 

2.7 

.2803 

2.7 

-.7171 

3.0 

.3366 

3.0 

-.71 37 

3.3 

.3960 

3.3 

-.6713 

3.6 

.4580 

3.6 

-.5896 

3.9 

.5221 

3.9 

-.4703 

4.2 

.5880 

4.2 

-.3171 

4.5 

.6551 

4.5 

-.1351 

4.8 

.7232 

4.8 

.0700 

5.1 

.7919 

5.1 

.2916 

5.4 

.8611 

5.4 

.5237 

5.7 

.9305 

5.7 

.7611 

6.0 

1.0000 

6.0 

1 . 0000 
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\ = 24.9 microinches 
M = 1.0 

6 = .006 inches 
1 = 6.0 inches 

With these, eq (3.26) becomes 

J 0 (2856 u f ) = 0 (3.38) 

With zero values of J 0 (x) presented in Table 7 and eq (3.38), a plot 
of fringe order as a function of in-plane displacement was determined' 
and presented on Figure 26. 

With this information, fringe order as a function of length along 
the beam for the first two modes of vibration can be determined. For 
the first mode of vibration, eq (3.36) and Figure 26 along with the 
fact 6= .006 inches and 1 = 6.0 inches was used and the results are 
plotted in Figure 27. For the second mode of vibration, eq (3.37) 
and Figure 26 was used and the results are plotted in Figure 28. Figure 
29 is a picture of the fringes for the first mode as they would appear 
in the transformed plane. Figure 30 is a picture of the fringes for the 
second mode of vibration. 

Point-by-Point Interpretation 
For the point-by-point interpretation, let 
\ = 24.9 microinches 
z = 40-.0 inches 
M * 1.0 

5 = .006 inches 


1 =6.0 inches 
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Therefore for n = 1 , eq (3.28) becomes 


.. _ 3.8 x 10" 4 

u f p 


(3.39) 


Thus, fringe spacing as a function of displacement was determined and 
presented in Figure 31. Plots of fringe spacing as a function of 
length along the beam for the first and second modes of vibration were 
also determined and are shown in Figures 32 and 33, respectively. 

Equipment for experimental verification of the theory presented 
in this section was unavailable. Thus, no experimental data is pre- 
sented for this section. 





3 
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Figure 26. Fringe Order as a Function of In-Plane Displacement 
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Figure 32. Fringe Spacing Along the Beam for the First Mode of Vibration 








IV. CONCLUSIONS 


The theory for the formation of fringes in the transform plane 
for double exposure and time averaging speckle interferometry was 
developed. Thus, in-plane displacement resulting from either static 
or dynamic loading can be determined. For the evaluation of the 
fringes, the whole field concept was used and the relationship be- 
tween fringes in the transform plane and displacement was derived. 

The technique is simple and needs no calibration and no laborious 
alignment of optical elements for the formation of the fringes. 

The sensitivity of the technique is dependent upon the resolving 
power of the film, the magnification factor and the location of the 
aperture in the transform plane. In addition, the magnitude of the 
displacement has to be larger than the speckle size which is dependent 
upon the aperture opening of the camera lens. 

In general, speckle interferometry has the ability to measure 
in-plane displacement for either the static or dynamic case within 
an accuracy of ^5 % for most of the displacement field. 
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I. INTRODUCTION 


1 .1 Basic Phenomena 

The relation between the elastic and electromaqnetic material 
properties was studied by Wertheim in 1844.( Annales de Chimie ,T. 

XII pp. 610-624). This interest was continued by Lord Kelvin who in 
1856 noticed that the resistance of wires changed as the wires were 
stretsfred [1]. This later led to the developement of the electri- 
cal-resistance strain gage. 

Since World War II strain gages have been used extensively by 
industry to determine strains in variously loaded members. The 
major disadvantage is the fact that three gages are required (rossette) 
to analyze a single point thus requiring many gages for a complete 
analysis [2]. 

Recently several optical stress analysis techniques have been 
developed tnat give a full field representation of displacement, 
strain or stress fields which utilize the wave nature of light as the 
mathematical model [3] . 

r - A cos (tit ++ } (1.1) 

where 

E - optic vector . tot - phase 

A - maximum amplitude $ - initial phase 

From Euler's identity 


exp(i9) = cosO + i sin0 


1 


( 1 . 2 ) 


2 


E becomes 

E* = A expEi (ut + <j>)3 

where it is understood that we are dealing with the real part of 
the exponential function. 

The intensity is 

I = pE-f* 

where E* is the complex conjugate , 'and- v is a constant. 

1 .2 Photoelasticity 

Photoel asticity uses the property of birefringence in certain 
materials [4]. When polarized light is incident on certain crystals 
it can be resolved into two components along the pricipal directions. 
One component is retarded with respect to the other and is referred 
to as the slow axis. The other is the fast axis. For a more rigor- 
ous discussion see [3] 3 or [5]. 

There are two basic arrangements for photoelasticity [5], the 
cross polariscope and the circular polariscope (Figure 1). 

In a crossed or standard polariscope, the two elements are a 
polarizer and an analyzer. Both are polarizing prisms that convert 
ordinary light into plane polarized light. There is no light output 
from the analyzer since it is '‘crossed" relative to the polarizer. 

If a birefringent model is placed between the polarizer and 
analyzer (Figure 2), then loaded, there is light output. 








SOURCE 
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The light from the polarizer can be described by 


K q exp(iwt) 


When light enters the model, it is resolved into two components 


Ep = ? cose XD - A 0 exploit) cosa 


xp 0 


xp 


E q = it sine x p = K q exp(.itot) sine xp 
When it leaves the model, the light vector becomes 


Ep i = sxp ( icut + 4 . ) cos e xp 
E i = K q expdrnt + <p ) sin 0 xp 


and when it leaves the analyzer 


^r “ E p' s,n0 xp * ^q 1 CQS0 xp 


which reduces to 


E r = A Q exp(iwt) sin 2s xp exp(i$ p ) - exp(i* q ) 


The intensity becomes 


I = yA 2 sin 2 2e xp [l - cosi <&] 


or 


I = u A 2 sin 2 2e xp sin 2 ( a <p ) 


(1.3) 


(K4) 


(1.5) 


( 1 . 6 ) 


(1.7) 


( 1 . 8 ) 


(1.9) 



V- \ 












* 


Figure 3. Circular Polar 
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t s = 1 / 2/2 A Q exp(iwc) (1.13) 

After leaving the quarter wave plate, the light vectors become 
Efi - /2/2 A q expi(wt + $.j?) 

? s , = /2/2 A q expi (ait + $ s ) (1.14) 

where 

$f = ir/2 + <£ s 

Upon entering the model, the light is resolved into two light vectors 
along the principal axes, p and q. They are 

Ep = E-, (cos(ir/4 - Q xp )) + ^" s i (sin(ir/4 - 0 ^)) 

f q “ ? f , (sin(,/4 - o xp J) - f s , (cos{ir/4 - 9 xp )} (1.15) 

and upon leaving the model, these light vectors become 

Epi = 1 A q /2 expi (wt + 4»s - ir/4 + 0 X p + 4> p ) 

~~T 

? q i = - A Q /2 expi (wt + $s - ir/4 + e xp + a ) (1.16) 

Note that 

exp ( i tt) = i 

2 

and that 

iexp(ie) = icose -sine . $ JJ 
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Upon entering the second quarter wave plate 


e s = V cos( i- v + K- sin( f ■ V 

gf = E pl Sint, - 0 xp ) - E q , cost - 8 xp ) 


(1.17) 


and after leaving the second quarter wave plate 

E c i = K n /2 expi (ait + 2$ s - t + 9 xp )[iexp(i$ n ) cos(^ 


o 

2 


4 


e xp> 


- exp (id) sin(^ - e )] 

9 4 X P 

If i = KJZ expi (tot + 2d - * + 0 Yn )[~exp(i<i> ) s1n(r - ? k J 
r s 4 X P P 4 


+ iexp(id )cos(ir - e )] 
9 4 xp 


0.18) 


Finally going into the analyzer 


C ■ t'/r- 1, i 


s 2 


(1.19) 


and after leaving the analyzer E r reduces to 

l* r = A Q expi (cot + 2* s + 2e xp )[exp(i<f. p ) - exp(i$ q )3 

The 1 ; intensity becomes 


( 1 . 20 ) 


1= ]i A 2(1- cos A <j) 
or 

I = uA2 sin bjJ 
0 2 

Fora fringe 


( 1 . 21 ) 


*10 


&*-n,,ir 0-22) 

T 11 

and from the stress optic law [5] , 

H. = f ff h(p - q) 

Where 

fcr is a stress-optic coefficient 
h is plate thickness 
(p - q) differences in principal stresses. 

So 

p - q « n.ji fc = 2 tmax (1.23) 

nr 

Equation 1.23 is the basic equation for photoelasticity. The 
information above, along with the isoclinic information, gives an 
accurate description of the shear stress field. 

Photoelasticity has many applications in stress analysis, for 
example it as: used as a tool in fracture mechanics»[7,8] for analyzing 
cracks. It can also be used in finding stress concentration factors 
for various structures and has application for wave propagation [9], 

1.3 Holography 

Holography takes advantage of the interference of two coherent 
light waves at a point on a film plane. Tne same mathematical model 
can be used as described in 1.1. Figure 4 shows the optical arrange- 
ment for holography. There are two beams, an object which illuminates 



LIGHT 

SOURCE 



Figure 4, Optical arrangement for holographic interferometry. 
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the object, and a reference beam which is directed to the film plane. 
This is the off-axis reference beam method developed by Leith and 
Upatnieks [10], 

Holographic interferometry is a double exposure technique. An 
exposure is made of the object before deformation, and an exposure 
is made after deformation. The surface displacement of a point is 
recorded as a phase variation, and fringes are produced. 

With the arrangement shown in Figure 4 [11] 


E r « A r expi (ait + e R ) 

t Q = A 0 expi (cut + e fl ) (1.24) 

The total li'ght vector for the first exposure is * 

0-2S) 

= ff R expi (wt + e R ) + A Q expi(wt + e Q ) (1.26) 

For the second exposure after deformation 

tg = A r expi (wt + e R ) + Aq expi (at + e Q + a 6 q ) (1.27) 

where A9 q is the change in the object wave due to deformation of 
the object by loading. 

The intensity at the film for the two exposures is 



£ * 



(1.28) 
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After developing the hologram, the film is then reconstructed 
(Figure 5), where 

^Recon = 1 A R exp ^ wt + 9 r) (1*29) 

From this, three distinct terms may be separated: a virtual term, 
a real image term, and an undiffracted term. 

The virtual term is used for analysis, and it reduces to p2j. 

^Virtual = 2 A R A o ^ * cosi9 ^ (1.30) 

As in photoelasticity , the&e term .can be interpreted physically. 

Ranson [12] shows (Figure 6) that theAQ term is related to a geometry 

change or 

&e = ZjH C (ps + po) • pp'] (1.31) 

X 

where 

ps is a unit vector from a point on the body to the light source, 

po is a unit vector from a point on the body to the film, 

pp' is the displacement vector. 

The final intensity expression becomes 

^Virtual = 2 A R A o cos (Mp s + P°) * PP 1 )] 0-32) 

X 

A fringe occurs when 1 = 0 or when cos&e = -1. 

With the use of a photopolymer, holography can be used to find 
natural frequencies in complex vibrating objects [13]. It is very 





Figure 5. Reconstruction of Hologram 
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Figure 6. Geometry of Fringe Formation 
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sensitive to out of plane motion, and can be used in wave propagation 
studies when a pulsed laser is used [11], 

Several optical techniques can be used to determine strains [12]. 
These are all based on geometry changes in a grid pattern either etched 
or projected onto the object. The most popular- of these is Moire' [15], 
Recently the speckle effect has received attention as a possible 
means of determining displacement components. This thesis will deal 
with this speckle effect in more detail. 


II. SHEARING SPECKLE INTERFROMETRY 
2.1 Static Loading Theory 

The optical configuration for shearing speckle interferometry 
is shown in Figure 7. Note the use of the wedge in the system. 
Because of this wedge, the focused image has superimposed on it a 
"shifted" image in the direction of the shea*". For a shear in the 
x-j direction, a point P is imaged at (x-j, Xg) and (x-j + Ax-j, Xg) in 

the film plane, and a point P-j is imaged at U-j, x g ) in the film 

plane. 

In order to obtain fringe data, a double exposure technique 
is used. A photograph is taken of the body in some reference position, 
The body is then deformed and another exposure is .superimposed on 
the same film. 

Ranson and Swinson [16] have shown that at the point (x^, Xg), 

the light amplitude for the first exposure will be given by 

E T i = E p + E p] 
where 


Ep~ A exp£i 9 p (x 1 ,Xg)] 

E p y A exp[i 9p-j(x^,Xg)] 

The total intensity for the first exposure is 


( 2 . 1 ) 


i- t * 

1 Ti Ti 


( 2 . 2 ) 
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i i v.~s? ( v •s-arw'l ; ■ ■—;* ! 



Figure 7, Optical configergtion for shearing speckle interferometry. 
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where denotes the complex conjugate 
or 

I-j - -EA exp[1 9 p (x r x 2 )l + I exp[i Sp^x^Xg)]} 

• {A exp[-ie p (x.j ,x 2 ) + A exp[-ie pl (x 1 a x 2 )]} 

= A 2 + A 2 exp[i8p(x^ ,x 2 ) - iSp.^XpXg}] + 

A 2 exp[i9 pl (x r x 2 ) - i0 p (x r x 2 )] + A 2 

= 2A 2 * 2ft 2 costSpCx-! »x 2 ) - 0 p1 (x 13 x 2 )] (2.5) 

For convenience we define Q = 0 p(x-j,x 2 ) - 0 p](x-j,x 2 ) 
so 

I-, = 2 ft 2 + 2ft 2 cose (2,6) 

Similarity for the second exposure 

Ej 2 = Epj + Ep-j i (2.7) 

where 

tp. = expi CSp -lie ] . 

' E pl , =ffexpl£9 pl -de pl J 

j A. © p is the change in phase due to deformation. The total intensity 

for the second exposure will be 


(2.3) 


(2.4) 
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I 2 = ^ T2 • f T2 * (2,9) 

= A" expi£e p + A.0p] + A expi£e pl + A.0 pl 3 • 

A exp(-i£e p + a 0 p 3 ) + I exp(-i[0 pl + A0 pl ]) (2.10) 

which reduces to 

I 2 = 2A 2 + 2A 2 cos (6 +A0) (2,11) 

where a 9 = A 0„ - a e„, 

P Pi 

So the total intensity expression for both exposures is 

t T = + I 2 (2,12) 

= 4A 2 + 2A 2 cose + 2A 2 cos(e +A.6) (2.13) 

Fringe data is obtained by taking the optical transform of the 
film as illustrated in Figure 8. For a transparency, the amplitude 
transmission function is linear for the ranges of interest or 
[16], 

gOO = a+bI T 


Where g(x) is the transmission function. The Fourier transform of 

gOO ^ 


[2. 14) 



s' 


Figure 8. Optical Configuration for a Fourier Transform. 
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0{w) = exp[-izw /2k] .Tg(x) exp(-iwx)dx (2.15) 

Where G(w) is the Fourier transform function of g(x) which is the 

light amplitude in the transform plane,/ ■ 

k = 2tt/x and 

w = kfi. . 
z 

Now define 

(2.16) 

(2.17) 

(2.18) 

G(w) = exp[~izw 2 ]/°{a + b[4A 2 + 2A 2 cos0 + 2X 2 cos(e -fee)] 

2k 

exp[-iwx]>dx (2.19) 

= exp£-izw 2 ]/*{a + b[4A 2 + 2A 2 cose(l + cosi.0) 

2k 

- 2A sine sime]} • exp[-iwx]dx (2.20) 

= exp[~izw 2 ][(a + 4bA 2 )<s‘w) + 2b4?C(w)(l +cosi9) 

2k 

- 2bA 2 S(w)sin a e * (2.21) 

The displacement information is related to the minimum of | g(w}|. 
or when 0 + cosie) = Q and sime = Q or 


<5(w) “ /^expC-iwxldx 
C(w) =* jocose exp[-iwx]dx 
S(w) = j“sine exp[-iwx]dx 
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^0 = (2n - 1 ) 7T n= 1 ,2,3,.. . 

As in other optical techniques, this 4.0 term can tie physically 
interpreted, Ranson and Swinson [16] have shown 
that for a wedge, the wedge angle a is related to the deviation 
angle 6 by (Figure 9) 

6 - (m ** 1 ) a ' (2.22) 

Where m is the index of refraction and the amount of shift or deviation 
M is 

M = On - T)®s 1 (2,23) 

Where s 1 is the distance from the lens to the focal point (Figure TO). 

When an object is illuminated by a single beam (Figure 11), the 
wedge causes points P and to be imaged at the same film location, 

The separation on the film plane isM* 

Consider now the phase variation of the two points P and P^, 

As in holography the displacement sectors are 

-3- I 

PP = u 1 e i 

P-jP] = u.e. - u.tx i +Ax.)e. (2,24) 

The phase change is now 

i0 2 = MP] S + P^] P-j 
A 




(2.25) 


A 9-. = 2tt C ps + po] • pp‘ 

1 X 

or for parallel unit vectors 

A9 = 2tt_[ps + po] * Cp-rP-j ' - PP 1 ) 

X 

and fringes occur when A e = Tr,3ir,,. ,(2n - 1 )tt 
so [16] 

(2n - l)ir = 2 tt (.1 + 1 n ) (§u, ax + au, Ay + au az) + 
x s 0 3X 9y az 

(m + ra ) (a_v ax + av. Ay + a_v az) + 

5 0 ax ay 3 z 

(n + n ) (3_w^ ax + 3_w_ Ay * a_w_ az) + 
s 0 ax ay az 

An experimental example of this technique is shown in section 


(2.26) 


(2.27) 

4.1. 
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2.2 Time Average 

Consider again the arrangement of Figure 7 . This time the 
object is caused to vibrate toward the camera, usually by a frequency 
generator or oscillator. A film record is made of the object during 
vibration. During the vibration, the object reaches two positions 
of zero motion. This effect is in essence the double exposure or the 
recording at two different configurations. The film record is pro- 
cessed and placed in a Fourier filter system as before. 

Again 




where 

E p = A((x 1 + x 1 , )>(x 2 + x 2 '})expi[e p (x r x 2 ) +*6^ 
and 


(2.28) 


t p i « A((x 1 + x 1 '),(x 2 + x 2 '))expi[ 0 pl (x 19 x 2 ):t Ae pl ] (2.29) 


Note that the amplitude A is a function of the frequency of vibration 
andA.9 p andAQpi are the relative phase changes due to vibration. 

The total intensity I is 


I = 'S<W • V* V )dt 

■ expi[e p +48^1 + A* expi[8 pl +I8 pl ] . 


(2.30) 
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A exp[-i(0 p +.A 0 p )] + A exp[-i(e pl + A9 pl )] 

= /J 2A 2 + Fexpi[ 9 p + A 0 p ** e pl -A0 pl ] + 

F exp[-n(e p +Ae p - 0 p -j -A0 pl )] 

* /g 2A 2 + A" 2 e<p[i(0 +A0)] + A 2 exp[-i(e -JA9 )] 

where 

Q = Q -e nl 
P pi 

A 0 = A0„ - A0„, 

P PI 

Recall g(x) = a+ bl 
which becomes 

g(x) = a + b/^ 2A 2 + A 2 exp[i(9 + A 0 )]+ A 2 exp[-i (0 + A8)]dt 
= a + b[2A* 2 t + A 2 exp(i 0 )/pexp(il 0 )dt + A 2 exp(-i0)/^exp 


As before 

G{w) = exp(jzv£)/g(x) exp(-iwx)dx 
2k 


t 

= exp(izw 2 )^ + bl^A^t + A^expCie)/- exp(iA8)dt 
2k 0 

f 

+ A 2 exp(“i0)/^ exp(-iA0)dt] } exp(-ixw)dx 

= exp(izw 2 )a <s(w) + 2A 2 bt <s(w) exp(izw 2 ) 

2k 2k 


+ bexpdswiJC^exptieJ/^expfii e)dt 


(2.31) 

(2.32) 

(2.33) 

(2.34) 

(2.35) 

- i Ca 0 )3 dt 

(2.36) 

(2.37) 

(2.38) 
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+ A exp(-l0)/^ exp(-iA ©)dt exp[-iwx]dx (2.39) 

As before the A 0 term can be expressed as 

A 0 = 2jt. (1 + l n )(5u A X + 3U Ay + 3U A z) + 

X ° 3 X ay 32 

(m * nu)(3wA x + 3w A y + 3 _wa z) + 
s u 3X 3y 3 z 

(n + n )(3WA x + 3 wa y + 3 wa z) 
s 3X 3y 3 z 

Now consider a specific optical arrangement (Figure 12). 

1 Q - 1 l s = cos e 

m Q = 0 m s = 0 

n Q = 0 n s = sin 8 

A x = Ay = 0 
So a 0 now becomes 

A 0 = 2TT, (1 + cos 0)<yj_ + singsw A z (2.40) 

X 32 32 

Assume 3w is negligible and u = u (z) sin t , where u is the 

3 Z 

maximum amplitude of vibration and to is the frequency of vibration. 

3u_ = 3 Uq sin tot (2.42) 

3 Z 3 Z 


For convenience define 
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B + exp(izw 2 ) s(w)[a + 2A 2 bt] (2.43) 

2k 

Now substituting Equations 2.40 - 2.43 into equation 2.39 


G(w) = B + bexp( izw 2 ) /[A 2 exp(i0)/g expi[2ll(1 + cosg)A z sun sinoitldt 
2k 4 3z 

+ A 2 exp(-ie)/pexp(~iC2J(l + cosfl)A z3uft simot1)dt exp(-iwx)dx 

3z (2.44) 

Recall that 


2ttJ o (x) = /^expdx sin©)de (2.45) 

if x. ~ 2ir (1 + cosg)3ugAz 

^ 3 Z 

0 = (lit 

d0 = lodt 
t-| = 2 tt / to 

then 


G(w) = B + bexp( izw 2 ) /tA 2 exp( i 0) J (2ir(l + cossJazsll) 

2k °~ n° 

+ ^expCio) 0 (2ir.(l + cosg)iz3_u)] exp(-iwx)dx (2.46) 

A 3Z 

Relative minimums occur when J Q (x) = 0 or when x = where ^ values 
are tabulated in Table 1. 

So 

2ir •(! + cosg)i z 3 u„ = 6 . 

T it 1 


(2.47) 
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TABLE 1 

ZERO VALUES FOR J Q (x) 


= 11.79 
= 14.93 

- 18.07 

- 21.21 
= 24.34 
= 27.49 
= 30.63 
= 33.78 
= 36.92 
* 40,06 


= 43.20 


i 


III. SINGLE BEAM SPECKLE INTERFEROMETRY 


3,1 Static Loading 

An object is illuminated by the laser light before deformation, 
and a photograph is taken with a conventional camera system as illust- 
rated in Figure 13. The object is then deformed, and again another 
photograph is taken superimposed on the first. 

For the first exposure, the light vector can be expressed 
as 


- A(x-j ,x 2 ) axpCie^x-j »x 2 )] (3,1) 

where 

(x-jjXg) - film plane coordinates 
sC.X'jaXg) - phase at film 
A(x-j,x 2 ) - amplitude factor 
The intensity expression for the first exposure is 

I = £, * E-j* (3.2) 

Now for the second exposure, again 

E 2 (x'-j,x‘ 2 ) =A'(x' 1 ,x , 2 ) exp[i 0 (x'-| ,x' 2 )] (3,3) 

and a similar expression exists for I 2 

Wilson [17] has shown that the total intensity at the film plane 
is 
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I T = A2(x r x 2 ) + A2(x‘ r x' 2 ) (3,4) 

Wilson further shows the relation between the intensity distribution 
to the displacement of the surface. 

Iy = A2(x-j »X 2 ) t A2(x, + U^sXg + U 2< j.) (3,5) 

where 

u-j f = Hu-j 

u 2f = Mu 2 

M = Magnification factor 
u-jjU 2 - displacements in x-j »x 2 directions 
The amplitude transmission function g(x-j,x 2 ) for a photographic 
film can be approximated as a linear function of intensity for the 
ranges of interest, or 

g(x=| ,x 2 ) = a f bl-j- (3,6) 

where a' and b are film constants. 

So 

gCx-j s x 2 ) = a + b£A2(x-j ,x 2 ) + A2(x-j + u-j^,x 2 + u 2 ^,)J 

Again the Fourier transform Gf the transmission function must 
be obtained. For convenience the x 2 coordinate will be suppressed or 

g(x) = a + b[A2(x) + A^(x + u^)] 


(3,7) 
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Therefore as in equation 2.15 

G(w) = exp(izw 2 ) /g(x) exp(~iwx)dx (3,8) 

2k 

Wilson has shown that G(w) reduces to 

G(w) = b exp(izw 2 ) F[A 2 ][1+ exp(-iwu^)] (3.9) 

2k T 

where 

> 

F[A 2 ] = /A 2 (x) exp(-iwx)dx (3,10) 

The intensity in the transform plane is expressed by 

I p = G(w) • G*(w) (3.11) 

So Ip finally reduces to 

Ip = 2b* F[A 2 ](1 + cos (wu f ) ) (3.12) 

and fringes are defined when 

1 + cos(wu^) = 0 (3,13) 

or wUf = (2n - l)ir (3.14) 

or u f * (n - i) \z (3.15) 

T 2 P 


An experimental example of this technique is shown in section 4.3. 
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3.2 Time Average 

Again, refer to Figure 13 for this arrangement. A frequency 
oscillator is used to induce vibration, a single beam of laser light 
is used to illuminate the object, and a camera records the motion. 

The Tight amplitude can be expressed as [17] 

f(x-j,x 2 ) = A(x-j + x'-pXg + x’ 2 ) exp[i(e(x-j,x 2 ) + Ae] (3,16) 

where 

x-,x 2 - film constants 

x 1 -I ,x' 2 - function of amplitude of vibration, 
frequency, time, and location 

0 - phase angle 

A ©- change in phase due to vibration 

The film exposure time t is' 

tI T - /J f(x 1 ,x 2 ) • E*(x r x 2 )dt (3.17) 

Again, suppress the coordinate notation so 

tlj = E(x) * E*(x)dt s /J A(x +x')dt (3.18) 

Physically x 1 is related to the frequency of vibration and the time or 

x' = u f (x 1 ,x 2 ) sin wt (3,19) 

where u^ is the maximum displacement on the film plane. 


Recal 1 


g(x) « a + bly 

= a + b, A 2 {x + Uj. sinwt)dt (3.20) 

t 0 T 

This amplitude transmission function is then Fourier transformed 


or 


G{w) = exp(izw 2 )/g(x) exp(-iwx)dx 
2k 

As before, this expression is evaluated and is a function of 
the Bessel Function as Wilson has shown [17], 


(3.21) 


G(w) = b_ exp( izwj) FtA 2 ] J (u f w) 
tw 2k 0 T 

where F is the Fourier Transform function. Zeros occur when 


(3.22) 


J Q (u f w) = 0 


or 


J p( 2 *PU f ) » 0 

X Z 

Refer to Table I for zeros of the Bessel function, 
example is shown in section 4.4. 


(3.23) 

An experimental 


IV. EXPERIMENTAL VERIFICATION 


4.1 Double Exposure Shearing Speckle Interferometry 


Consider a cantilever beam with a shift in the z direction only 
(Figure 14). Equation 2.27 reduces to [16]. 



so 

(n-1/2 ) \ * [(1 + cose)H + sin ^|]Az (4.2) 


For a cantilever beam with tip deflection <s» then 


u ■ (x2 - tf)- + jiz 2 ] (4-3) 

w - %xy2 - (Iz - #)x] (4.4) 

13 c 

M- 3fi r -v(x^ - y2) - ~z* + Iz] (4.5) 

3Z "Ta S 2 

= %-lx fzx] (4.6) 

Sz 13 


Along the center of the beam y=0 and x=.125 in. 
Then 
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FRONT 


5 



A =24. 91 x 10 
65 =0.0899 in. 
o =0.37 
b =0.25 


Figure 14. Experimental arrangement for double exposure shearing 
speckle interferometry. 
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= (-.201 - 34.7z 2 + 416. 6z) x 10" 6 (4.7) 

= (-52.1 4- 8.68z) x 10" 6 (4.8) 

oZ 


By substituting 4.7 and 4.8, plus the other information, into 
equation 4.2, a theoretical fringe equation can be obtained as a function 
of beam length z. 


(n-l)x= [1.942(-.201 - 34.7z 2 + 416. 6z) x 10“ 6 


+ 0.335 (-52.1 + 8.68z) X 10" 6 ] (.0899) 

(4-9) 

or 


n = .436 + 2.93z - .244z 2 

(4.10) 


A shearing speckle photograph was taken and filtered, which 
is shown in Figure 15. 

Using Equation (4.10), a theoretical fringe plot was made and is 
shown in Figure 16. Using the shearing speckle photograph, the 
actual fringe order as a function of- distance can be found and is 
plotted on the theoretical curve of Figure 16. The agreement is 
obvious. 

4.2 Shearing Speckle Interferometry Vibration 

Again look at Figure 12 and let the object be a cantilever beam 
with the configeration shown in Figure 17. Nowacki [18] presents a 
solution for the amplitude of vibration along the beam as 

w p (z,) = cDJU.z) - v( u) (4.n) 

r TfeT 



I * ■ 



Figure 15. Shearing speckle photograph after Fourier filtering. 







FRINGE 0 


X 



DISTANCE (IN.) 


Figure 16. Comparison of theoretical curve and experimental points for double 
exposure shearing speckle interferometry. 
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ill! 


where 

W is the amplitude 
C is a constant 

u( i r z) = coshy; - CQS^z 
2 

V(iz) = sinhx 2 - sinX z 

r r 

2 

s( 3 r ) = 1. (coshB r + cose r ) 

T(P r ) - \ {sinhB r + sin& r ) 

is wavelength due to vibration 

e r = V 

The transcendental equation- for a vibrating Cantilever beam is .[19] 

cosh6 r cosa r = -1 (4.12) 

The solutions are tabulated in Table 2. For the first mode of 
vibration r=l and at 2=1 

«,(!) = 6 = C[ U(e 1 ) - f[|4v(8 1 )] (4.13) 

So 


W.J = U(X 1 2)- .734 V(X 1 2) 


(4.14) 


for the first mode. 
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TABLE 2 


Solutions for Transcendental 
Equation for Vibrating 
Cantilever Beam 


1.875 

4.694 

7.855 


10.996 
14. 137 
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Taking the derivative of W-j (1) the |~ term can be found. 


1.9 u „ Xi.[(sinhU,,z) + sinU,,z)) - 

5 9z ** 2 1 1 

l.734lcoshU*j,z) - cos(X 1# z))3 (4.15) 


A cantilever beam with' dimensions and properties shown in 
Figure 17 was vibrated in the first mode, with two different tip 
deflections, and shearing speckle photographs were taken. These 
photographs were Fourier filtered and drawings are shown in Figure 
18. 

From Equation (2.47) 


1 3 u _ $ . , 

Tal " 1 x 


(4.16) 


2ti(1+co5£)Az. 

From the theoretical equation (4.15) a curve is drawn in Figure 
19; the data points are plotted on this curve. The agreement is obvious. 



= 7 in. 

= .75 in 
£>z= 0.2456 in. 

50 cycles per second 


FREQUENCY 

OSCILLATOR 



Figure 17. Experimental arrangement for vibration shearing speckle. 


ze 9 



Figure 19. Plot of theoretical curve and experimental points for vibration 
shearing speckle. OS=1.2xlO" 3 ©fi=8.8xl(T 3 
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4.4 Single Beam Speckle Vibration 

Figure 23 shows the experimental configuration used. Again 
using the solution by Nowacki [18] for the first mode 


W «■ UU r z) - *734 VU v z) (4.20) 

5 

This gives a theoretical curve for the vibration amplitude. 

Recall Equation (3.23) 

Jq ( 2irpu,j>/ Az) = 0 (4.21) 

For the point-by-point data reduction (Figure 24), the distance 
between the first two fringes is the zero of the Bessel function 
giving 

2irpu^/AZ - 2.40 (4.22) 

For the distance between the second outer fringes, the second 
zero of the Bessel function gives 

2irpu^./ AZ = 5.52 (4.23) 

Table 3 gives the results of the data reduction. Notice at 
z = 4.95 in. through 6.67 in. , the second outer fringes were used, 
thus the change in the constant relating displacement to fringe 
spacing. 

Figure 25 shows a plot of the experimental and the theoretical 
displacement curve. Again the agreement is obvious . 
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Figure 24. Sample of Point by Point Data. 
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DISTANCE ALONG BEAM Z (IN.) 

Figure 25. Theoretical curve and experimental points for vibration single 
beam speckle. 
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Model 
Lacati on 

2.86 

3.14 

3.33 
3.62 
3.81 
4.00 
4.19 
4.38 
4.57 
4.76 
4,95 

5.14 

5.33 
5.52 
5.71 
5.90 
6.09 
6.29 


TABLE 3 

DATA REDUCTION FOR VIBRATION SPECKLE 

X Model u.p Model 


4.50 

1.56xl0~ 4 

3.56 

1.98xl0“ 4 

3.14 

2- 24x10 ~ 4 

2.95 

2.39xl0" 4 

2.78 

2.53xl0" 4 

2.57 

2.74xl0~ 4 

2.5 

2.82xl0 -4 

2.38 

2.96xl0" 4 

1.92 

3.67xl0~ 4 

1.85 

3. 80xl0~ 4 

4.08 

3.97X10" 4 

3.89 

4.16xl0~ 4 

3.60 

4.50X10' 4 

3.39 

4. 78x10 ~ 4 

3.20 

5.06xl0" 4 

3.10 

5.23xl0" 4 

2.90 

5.59xl0' 4 

2.76 

5.87X10** 4 


e.ooxio ” 4 


6.47 

6.67 


2.70 



V. CONCLUSIONS AND RECOMMENDATIONS 


The speckle effect has been noticed for several years by many 
investigators. The theory for shearing speckle and single beam 
interferometry was developed. 

Shearing speckle, like holography, is dependent upon a’ phase 
change. In general, for a shear in one direction at least two 
derivations of displacement terms are combined in the fringe pattern. 
Experimentally shearing speckle is simple to set up, and for double 
exposure the fringes are easy to transform. For time average ana- 
lysis, the fringes are difficult to transform. For this type of 
experimental work time average holography yields much clearer fringes, 
and with photopolymer recording, natural frequencies of complex 
vibrating objects can be found. 

Single beam speckle interferometry is dependent upon an ampli- 
tude change. This technique is very simple to set up experimentally, 
yielding two in-plane displacement components. For time average, 
single beam speckle can be used to determine the amplitude of motion, 
even for thin structures. 
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